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TAIL APPROXIMATIONS OF INTEGRALS OF GAUSSIAN 
RANDOM FIELDS 

By Jingchen Liu^ 
Columbia University 

This paper develops asymptotic approximations of P{Jj, e^'*' dt > 
6) as 6 — > oo for a homogeneous smooth Gaussian random field, /, 
living on a compact d-dimensional Jordan measurable set T. The in- 
tegral of an exponent of a Gaussian random field is an important 
random variable for many generic models in spatial point processes, 
portfolio risk analysis, asset pricing and so forth. 

The analysis technique consists of two steps: 1. evaluate the tail 
probability P(/_ e^''^ dt > b) over a small domain H depending on b, 
where mes(H) ^0 as b ^ oo and mes(-) is the Lebesgue measure; 
2. with H appropriately chosen, we show that P^j'^e^'^^^ dt > 6) = 
(1 + o(l)) mes(r) X mes"^(H)P(4 e-*^'*' dt > b). 

1. Introduction. We consider a Gaussian random field living on a d-di- 
mensional domain T C i?^, {f{t):t G T}. For every finite subset {ti,..., 
tn} C T, (/(ti), . . . , is a multivariate Gaussian random vector. The 
quantity of interest is 



T 

as 6 — oo. 

The motivations of the study of e-^'-*-' dt are from multiple sources. We 
will present a few of them. Consider a point process on associated with 
a Poisson random measure {Na}a&b with intensity \{t), where B repre- 
sents the Borel sets of R^. One important task in spatial modeling is to 
build in dependence structures. A popular strategy is to let f{t) =logA(t), 
which can take all values in i?, and model f{t) as a Gaussian random field. 
Then, e-''^*) dt = E{N{A)\\{-)) for all AeB. With the multivariate Gaus- 
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sian structure, it is easy to include linear predictors in the intensity process. 
For instance, [19] models f{t) = U{t) + W{t), where U{t) is the observed 
(deterministic) covariate process and W{t) is a stationary AR(1) process. 
Similar models can be found in [18, 23, 40] which are special cases of the 
Cox process [21, 22]. Such a modeling approach has been applied to many 
disciplines, a short list of which is as follows: astronomy, epidemiology, ge- 
ography, ecology, material science and so forth. 

In portfolio risk analysis, consider a portfolio consisting of equally weighted 
assets (5i, . . . , S'n). One stylized model is that (log^i, . . . ,log5„) is a mul- 
tivariate normal random vector (cf. [7, 12, 25, 26, 30]). The value of the 
portfolio S = Si is then the sum of correlated log-normal random vari- 
ables. If one can represent each asset price by the value of a Gaussian random 
field at one location in T, that is, log Si = f{ti). As the portfolio size tends 
to infinity and the asset prices become more correlated, the limit of the unit 
share price of the portfolio is liuin^ooS/n = JrpC^^^^ dt. For more general 
cases, such as unequally weighted portfolios, the integral is possibly with 
respect to some other measures instead of the Lebesgue measure. 

In option pricing, if we let S{t) be a geometric Brownian motion (cf. [14], 
Chapter 5 of [27], Chapter 3.2 of [29]), the payoff function of an Asian option 
(with expiration time T) is a function of Jq S{t) dt. For instance, the payoff 
of an Asian call option with strike price K is max (/o S{t)dt-K,0); the 
payoff of a digital Asian call option is I{Jf^ S{t) dt > K). 

We want to emphasize that the extreme behavior of jrpeJ^^^ dt connects 
closely to that of supj^ /(O- show in Theorem 1, with the thresh- 

old u appropriately chosen according to 6, the probabilities of events 
{jj,ef^^^ dt > b} and {suprpf{t)dt > u} have asymptotically the same de- 
caying rate. It suggests that these two events have substantial overlap with 
each other. Therefore, we will borrow the intuitions and existing results on 
the high excursion of the supremum of random fields for the analysis of 
/^e/W dt. 

There is a vast literature on the extremes of Gaussian random fields 
mostly focusing on the tail probabilities of supj^/(t) and its associated ge- 
ometry. The results contain general bounds on P(sup^/(t) > b) as well as 
sharp asymptotic approximations as 6 — )• oo . A partial literature contains 
[13, 16, 17, 31-33, 35, 37]. Several methods have been introduced to ob- 
tain asymptotic approximations, each of which imposes different regularity 
conditions on the random fields. A few examples are given as follows. The 
double sum method [34] requires expansions of the covariance function and 
locally stationary structure. The Euler-Poincare Characteristic of the ex- 
cursion set [x(^b)] approximation uses the fact that P(M > 6) E^xi^b)), 
which requires the random field to be at least twice differentiable [1, 5, 38]. 
The tube method [36] uses the Karhunen-Loeve expansion and imposes dif- 
ferentiability assumptions on the covariance function (fast decaying eigen- 
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values). The Rice method [9-11] represents the distribution of M (density 
function) in an imphcit form. Recently, the efficient simulation algorithms 
are explored by [2, 3]. These two papers provided computation schemes that 
run in polynomial time to compute the tail probabilities for all Holder con- 
tinuous Gaussian random fields and in constant time for twice differentiable 
and homogeneous fields. In addition, [4] studied the geometric properties of 
a high level excursion set for infinitely divisible non- Gaussian fields as well 
as the conditional distributions of such properties given the high excursion. 

The distribution of / e-^*-*^ dt for the special case that fit) is a Wiener pro- 
cess has been studied by [27, 39]. For other general functionals of Gaussian 
processes and multivariate Gaussian random vectors, the tail approximation 
of the finite sum of correlated log-normal random variables has been studied 
by [8]. The corresponding simulation is studied in [15]. The gap between 
the finite sums of log-normal r.v.'s and the integral of continuous fields is 
substantial in the aspects of both generality and techniques. 

The basic strategy of the analysis consists of two steps. The first step 
is to partition the domain T into n small squares of equal size denoted 
by Ai, i = 1, . . . ,n, and develop asymptotic approximations for each pi = 
P(X4.e-^W dt>b). The size of Ai will be chosen carefully such that it is valid 
to use Taylor's expansion on f{t) to develop the asymptotic approximations 
of Pi. The second step is to show that P{Jj,e^^^^ dt>b) = (1 + o{l)) XliLi^'*- 
This implies that when computing -P(/y e^^^'^ dt > b), we can pretend that all 
the J^ e^^^^dfs are independent, though they are truly highly dependent. 
The sizes of the Aj's need to be chosen carefully. If Ai is too large, Taylor's 
expansion may not be accurate; if Ai is too small, the dependence of the 
fields in diff'erent Aj's will be high and the second step approximation may 
not be true. Since the first step of the analysis requires Taylor's expansion 
of the field, we will need to impose certain conditions on the field, which 
will be given in Section 2. 

This paper is organized as follows. In Section 2 we provide necessary 
background and the technical conditions on the Gaussian random field in 
context. The main theorem and its connection to asymptotic approximation 
of P{supj' f{t) > b) are presented in Section 3. In addition, two important 
steps of the proof are given in the same section, which lay out the proof 
strategy. Sections 4 and 5 give the proofs of the two steps presented in 
Section 3. Detailed lemmas and their proofs are given in the Appendix. 

2. Some useful existing results. 

2.1. Preliminaries and technical conditions for Gaussian random field. 
Consider a homogeneous Gaussian random field, fit), living on a domain T. 
Denote the covariance function by 

C(i-^) = Cov(/(s),/(t)). 
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Throughout this paper, we assume that the random field satisfies the fol- 
lowing conditions: 

(CI) / is homogeneous with Ef{t) = and Ep{t) = 1. 

(C2) / is almost surely at least three times continuously differentiable 
with respect to t. 

(C3) T is a d-dimension Jordan measurable compact subset of R'^. 

(C4) The Hessian matrix of C{t) at the origin is — /, where / is a d x d 
identity matrix. 

Condition (CI) imposes unit variance. We will later study Jj,e'^^^^^ dt and 
treat a as an extra parameter. Condition (C2) implies that C{t) is at least 6 
times differentiable. In addition, the first, third and fifth derivatives of C{t) 
evaluated at the origin are zero. For any /(t) such that AC(0) = S and 
|S| > 0, (C4) can always be achieved by an affine transformation on the 
domain T by letting f{t) = f{T}/h) and 

/ e-/(*) dt= [ e^f'-^"'"'^ dt = / e'^^(^) ds, 

Jt Jt J{s:Y:-^/^seT} 

where for a symmetric matrix S we let S^/^ be a symmetric matrix such 
that sV2sV2 = s. 

For a > 0, let 

(1) Ia{A)= [ e'^^Wdt 

J A 

for the Jordan measurable set AcT. Of interest is 

P{lAT)>b) 

as 5 — 7- oo. Equivalently, we may consider that the variance of / is cr^. How- 
ever, it is notionally simpler to focus on a unit variance field and treat a as 
a scale parameter. 

We adopt the following notation. Let "(9" and "A" denote the gradient 
and Hessian matrix with respect to t, and "(9^" denote the vector of second 
derivatives with respect to t. The difference between "A" and "(9^" is that, 
for a specific t, A/(t) is a d x d symmetric matrix whose upper triangle 
entries are the elements of d"^ f{t) which is a ((i((i-|- l)/2)-dimensional vector. 
Let dj denote the partial derivative with respect to the jth component of 
t = (ti, . . . ,td). We use similar notation for higher order derivatives. For b 
large enough, let u be the unique solution to 

(2^/a)'^/\-'^/2g<x«^^_ 

The uniqueness of u is immediate by noting that the left-hand side is mono- 
tone increasing with u for all u > d/{2a). In addition, we use the following 
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notation and changes of variables: 

^^l{t) = -i^lCit),...,^dCit)), 

fi2{t) = {dlC{t),i = 1, . . . , d; i = 1, . . . , d, j = i + 1, . . . , d), 

^02 = ^20 = ^2(0), /(0) = u-u;, df{0) = y, 
d^fiO) = ufio2 + z, Af{0) = -ul + z. 

The vector /U20 contains the spectral moments of order two. Similar to A/(0) 
and 9^/(0), z is a symmetric matrix whose entries consist of elements in z. 
We create different notation because we will treat the second derivative 
of / as a matrix when doing Taylor's expansion and as a vector when do- 
ing integration. As stated in condition (C4), we have AC(0) = — /. Equiv- 
alently, 5/(0) is a vector of independent unit variance Gaussian r.v.'s. We 
plan to show that in order to have e-^^*) dt > b, supj^ f{t) needs to reach 
a level around u. The distance between /(O) and u is denoted by w. In 
addition, since (/(O), 5^/(0)) is jointly independent of 5/(0), the distribu- 
tion of (9/(0) is unaffected even if /(O) reaches a high level. Further, the 
covariance between /(O) and 9^/(0) is /i20- Given /(O) = u, the conditional 
expectation of 9^/(0) is it/xo2- The distance between 9^/(0) and this condi- 
tional expectation is denoted by vector z. 

A well-known result (see, e.g.. Chapter 5.5 in [5]) is that the joint distri- 
bution of (/(O), 0^/(0), 9/(0), /(t)) is multivariate normal with mean zero 
and variance 

/I m c{t) \^ 

M02 Ai22 l^lj{t) 

0/ fij{t) 

\C{t) ^2(0 W(i) 1 / 

where fJ.i{t), H2{t) and 1^20 = fJ^Q2 defined previously. The matrix 1^122 is 
a d{d + l)/2 by d{d + l)/2 positive definite matrix and contains the 4th- 
order spectral moments arranged in an appropriate order. Conditional on 
/(O) = u — w, 5/(0) =y and A/(0) = —ul + z, f{t) is a continuous Gaussian 
random field with conditional expectation 



(2) E{t) = iu-w,un20 + z\y') 








where 

(3) 



1 ^20 
^02 ^22 

Note that Ufi20 + ^ is the vector version of —ul + z. Therefore, conditional 
on (/(0),9/(0)~'', (9^/(0)^) = {u — w,y~^ ,ufi2o + z~^), we have representation 

f{t) = E{t)+git), 
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where g{t) is a Gaussian random field with mean zero and 

E{t) = E{f{t)\f{0) = u-w, df{0) = y, 02/(0) = n^02 + z), 

whose form is given in (2). Since C{t) is six times differ entiable, E(t) is at 
least four times differ entiable. Using the form of E{t) in (2) and T in (3), 
after some tedious calculations, we have 

E{0) =u-w, dE{0) = y, 

(4) AE{0) = -uI + z, dfjkE{0)=y^d^jkfim, 

In order to obtain the above identities, we need the following facts. The first, 
third and fifth derivatives of C{t) evaluated at are all zero. The first and 
second derivatives of C(t) are contained in iJ,i{t) and /i2(i)- We also need to 
use the fact that 

_ 1 - M20^22Vo2 1 - A^20^22Vo2 

-1 -1 -1 ■ 

^22 /^02 -1/^22/^02/^20/^ 

— 1 /^22 ~'~ —1 

\ 1 — /^20/^22 /^02 1 — /i20/^22 /^02 / 

With the derivatives of E{t), we can write 

(5) E{t) = u-w + y^t + \t^{-ul + z)t + gsit) + g^it) + R{t). 
If we let t = {ti,... then 

(6) 93 (t) = ^ j;4,i?(0)tit,t,., 94 (t) = ^ E dfjkiE{0)Utjtkti, 

i,j,k ij,k,l 

and R{t) is the remainder term of the Taylor expansion. The Taylor expan- 
sion of E{t) is the same as f{t) for the first two terms because g{t) is of 
order 0(|tp). It is not hard to check that Var(g(t)) < c\t\^ for some c > 
and \t\ small enough. 

2.2. Some related existing results. For the comparison with the high ex- 
cursion of supy f{t), we cite one result for homogeneous random fields, which 
has been proved in more general settings in many different ways. See, for 
instance, [5, 9, 34]. This result is also useful for the proof of Theorem 1. 
For comparison purpose, we only present the result for the random fields 
discussed in this paper. 

Proposition 1. Suppose Gaussian random field f satisfies condi- 
tion (C1)-(C4). There exists a constant G such that 

P(supf{t)>u) = (l + o(l))Gmes(r)n'^P(/(0) >u) 

as u — )• oo . 
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We also present one existing result on the tail probability approxima- 
tion of the sum of correlated log-normal random variables which provides 
intuitions on the analysis of /^C'^^*^ dt. 

Proposition 2. Let X = {Xi, . . . be a multivariate Gaussian ran- 
dom variable with mean /i and covariance matrix S, with det(S) > 0. Then, 

(n \ n 

^e^'>6 =(l + o(l))^P(e^'>6) 
1=1 / i=i 

as 6 — oo . 



The proof of this proposition can be found in [8, 28]. This result implies 
that the large value of Y17=i is largely caused by one of the XiS being 
large. In the case that XiS are independent, Proposition 2 is a simple corol- 
lary of the subexponentiality of log- normal distribution. Though the Xj's are 
correlated, asymptotically they are tail-independent. The result presented in 
the next section can be viewed as a natural generalization of Proposition 2. 
Nevertheless, the techniques are quite different from the following aspects. 
First, Proposition 2 requires E to be nondegenerated. For the continuous 
random fields, this is usually not true. As shown in the analysis, we indeed 
need to study the sum of random variables whose correlation converges to 1 
when b tends to infinity. Second, the approximation in (7) is for a sum of 
a fixed number of random variables. The analysis of the continuous field 
usually needs to handle the situation that the number of random variables 
in a sum grows to infinity as 6 — )■ oo. Last but not least, to obtain approxima- 
tions for e-^W dt>b), one usually needs to first obtain approximations 
for P(/^ e^(*) dt > b) for some small domain C T. We will address all 
these issues in later sections. 

For notation convenience, we write = 0(6„) if there exists a con- 
stant c > independent of everything such that < cbu for all n > 1 , and 
au = o{bu) if CLu/bu —^0 as u — )■ oo and the convergence is uniform in other 
quantities. We write = @{bu) if = 0{by) and b^ = 0{au)- In addition, 
we write X^ = Op(l) if Xu A as u — )■ oo. 

3. Main result. The main theorem of this paper is stated as follows. 

Theorem 1 . Let f be a Gaussian random field living on T C R*^ satis- 
fying (C1)-(C4). Given a > 0, for b large enough, u is the unique solution 
to equation 

(8) (^Y\-d/2^^u^^_ 
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Then, 

p(^j e^-^(*) dt>h^ = (l + o(l))//mes(T)u'^-iexp(-nV2) 

as 6—7-00, where mes(T) is the Lehesgue measure of T , 

|r|-V2 det(/i22)^/2e((i/8)i^^22i+{i/8)E.9f.«C{0))/<x2 

^= (27r)('^+l)('^+2)/4 

(9) 

y_Rd(d+l)/2 \ 2 



X 



1 — A*20/^22^A*02 



r is defined in (3), fi20, /^02; ^22 ^'"e defined in the previous section and 

d d(d-l)/2 

Remark 1. The integral in (9) is clearly in an analytic form. We write 
it as an integral because it arises naturally from the derivation. 

Corollary 1. Let f be a Gaussian random field living onT C R'^ satis- 
fying (C1)-(C4). Adopting all the notation in Theorem 1, let b = b{27r / a)~^^'^ 
and 

(10) u = + A log f i^") + ^ ' ^^/^^ 



a 2a \ o- J \2 J cr log 6 
Then, 

P{Ia{T) > 6) = (1 + o{l))Hm^es{T)u'^-^ e^^{-u^/2). 

Proof. The result is immediate by the Taylor expansion on the left- 
hand side of equation (8) and note that u — u = o{u~^). □ 

As we see, the asymptotic tail decaying rates of sup^^ f{t) and J^e^f^^Ut 
take a very similar form. More precisely, 



p(^j e'^-^W dt>b^= e(l)P(^sup/(t) > 



with u and b connected via (8). This fact suggests the following intuition 
on the tail probability of I„{T). First, the event {/^(T) > 6} has substan- 
tial overlap with event {sup2^/(t) > u}. It has been shown by many studies 
mentioned before that given u sufficiently large {swprp f{t) > u} is mostly 
caused by just a single f{t*) being large for some t* G T. Put these two 
facts together, {lo-(T) > 6} is mostly caused by {f{t*) > n}, for some t* £T 
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not too close to the boundary of T. Therefore, conditional on {/^(T) > b}, 
the distribution of f{t) is very similar to the distribution conditional on 
{supji/(t) > u}. Of course, these two conditional distributions are not com- 
pletely identical. The difference will be discussed momentarily. Now we per- 
form some informal calculation to illustrate the shape of f{t) given f{t*) = u. 
Thanks to homogeneity, it is sufficient to study t* = 0. Given /(O) = u, 



Since C{t) is 6 times differentiable, dC{0) = and AC(0) = — /, we obtain 
E{f{t)\f[0) =u) ^ u — u- t^t/2. For the exact Slepian model of the random 
field given that / achieves a local maximum at t* of level u, see [6]. Note 
that for b large. 



In Theorem 1, this is exactly how u is defined. As shown in Figure 1, 
the three curves are exp{E{f {t)\f (t*) = u)} for different t*'s. Given that 
{supji/(t) > u}, these three curves are equally likely to occur. 

Second, as mentioned before, the conditional distributions of f(t) are dif- 
ferent given {If^{T) > b} or {supj^ f{t) > u}. This is why the two constants 
in Theorem 1 (H) and Proposition 1 (G) are different. The difference is 
due to the fact that the symmetric difference between {sup^ f{t) > u} and 



E{f{t)\f{0) = u)=uC{t). 




is approximately equivalent to 




t=0 



t 



Fig. 1. One- dimensional example. 
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{Jj,e^^^^ dt > b} is substantial though their overlap is significant too. Con- 
sider the following situation that contributes to the difference. sup'p/(t) is 
slightly less than u [e.g., by a magnitude of 0{u~^)]. For this case, Ia{T) 
still has a large chance to be greater than b. For this sake we will need 
to consider the contribution of A/(0). As is shown in the technical proof, if 
t* = argsup/(t) = and <9/(0) = 0, then a sufficient and necessary condition 
for Ia{T) > 6 is that 

/(O) + ^ Tr(n-^A/(0) + I)>u + o{u-^), 
la 

where TV denotes the trace of a squared matrix. Note that conditional on 
/(O) = n, ^(A/(0)|/(0) =u) = -ul. Therefore, A/(0) + ul is of size 0(1). 
One well-known result is that the trace of a symmetric matrix is the sum of 
its eigenvalues. Let Aj be the eigenvalues of A/(0) + / = z/u. Then, the 
sufficient and necessary condition is translated to /(O) + ^ Z^^Li \ > u. This 
also suggests that, conditional on I(j{T) > b, w = /(O) — ti is of size 0{u~^). 
This forms the intuition behind the proof of Theorem 2. 

The proof of Theorem 1 consists of two steps presented in Sections 3.1 
and 3.2, respectively. Each of the two steps is summarized as one theorem. 

3.1. Step 1. Construct a cover of T, {^1, A„}, such that T C |jr=i 
Each Ai is a closed square, mes{Ai n Aj) = for i ^ j. Because T is Jor- 
dan measurable, as supjmes(Aj) — )■ 0, nies(U"^;^ Ai) — mes(T) — )• 0. To sim- 
plify the analysis, we make each Ai of identical shape and let Ai = {ti + s : 
s G [0,0;,]'^}. The size of the partition n and choice of depend on the 
threshold b. The first step analysis involves computing the integral pi = 
P{J^. e"^'^^^ dt > b). Because / is homogeneous, it is sufficient to study pi. 

The basic strategy to approximate pi is as follows. Because / is at least 
three times differentiable, the first and second derivatives are almost surely 
well defined. Without loss of generality, we assume that € ^i. Conditional 

on (/(0),a/(0),A/(0)),/(i) = /(0) + a/(0)Tt+ltTA/(0)t + 53(t) + 54(t) + 
R{t) + g{t), where g{t) is a Gaussian field with mean zero and variance of 
order 0(|t|^). Then, 

W dt>b^ 

■) 

^„[u-w+y'Tt+il/2)t'T{-uI+7.)t+g3{t)+gi(t)+R{t)+g{t)] 

li 

• b I dwdydz, 
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where h{w,y,z) is the density function of (/(O), 3/(0), A/(0)) evaluated at 
(u — w,y,ufi20 + -2), which is a multivariate Gaussian random vector. Let u 
be defined in (8). There exists a S > (small enough) such that if we let Ai's 
be squares of size e = 0{l)u~^^'^~^^ and, hence, n = 0(1) mes(r)ii'^/^~'^^, the 
asymptotics of pi can be derived by repeatedly using Taylor's expansion and 
evaluating the integral on the right-hand side of (11). The main result of 
this step is presented as follows. It establishes a similar result to that of 
Theorem 1 but within a much smaller domain. 

Theorem 2. Let f be a Gaussian random field living in T satisfying 
conditions (C1)-(C4). Let Ai = 'E^ = {t: |t|oo < e}, where \t\oo = maxj|tj|. 
Let u and H be defined in Theorem 1. Without loss of generality, assume 
"Bi; cT with e = ku^~^^'^ for some 5 small enough and u large enough. Then, 
for any k > 

p^=p{E,) = p(^j^ e^f^'Utyb^ =(l + o(l))i/mes(H,)n"'-^e-"'/2 
as 6 — oo . 

The proof of this theorem is in Section 4. We will then choose each to 
be of the same shape as H^. Then, all the pi's are identical. 

3.2. Step 2. The second step is to show that with the particular choice 
of in the first step, P(/j,e'^^W dt>b) = (l + o(l)) Yji=iPi- We first present 
the main result of the second step. 

Theorem 3. Let f be a Gaussian random field satisfying conditions (Cl)- 
(C4) and e be chosen in Theorem 2. Let k G Z"^ and k = 2ke-|-H£. Further, 
let C" = {k : Hg^k C T} and C+ = {k : H^^k n T / 0}, then 



and 



p(lJ u H,,k) >&) = (i+o(i)) E p(u^e,^: 
u s.,k) = (i+o(i)) E nu^sM) 



>b) 



>b). 



^kec- ' ' kec- 

We consider 



Vec+ ' kec+ 



as a sum of finitely many dependently and identically distributed random 
variables. The conclusion of the above theorem implies that the tail distribu- 
tion of the sum of these dependent variables exhibits the so-called "one big 
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jump" feature — the high excursion of the sum is mainly caused by just one 
component being large. This result is similar to that of the sum of correlated 
log-normal r.v.'s. Nevertheless, the gap between the analyses of finite sum 
and integral is substantial because the correlation between fields in adjacent 
squares tends to 1. For finite sums, the correlation is always bounded away 
from 1. The key step in the proof of Theorem 3 is that the e defined in 
Theorem 2, though tends to zero as 6—7-00, is large enough such that the 
one-big-jump principle still applies. We will connect the event of high excur- 
sion of -^o-(UkeC+ '^£,k) to the high excursion of supy^^^^ =^ ^ f{t) and apply 
existing results on the bound on the supremum of Gaussian random fields. 
A short list of recent related literature on the "one-big-jump" principle and 
multivariate Gaussian random variables is [8, 24, 28]. 

With the preparation of the two steps, we are ready to present the proof 
for Theorem 1. 



Proof of Theorem 1. Prom Theorem 2, 

^ P{UE,^^)>b) = {l + oil))Hmes( (j H,,k) ^'^"'e-^Vs 
kec+ \ec+ ' 

P{Iai^e,k)>b) = {l + oil))Hmes( [j H,,k] u^-^e-^'/^ 
Therefore, thanks to Theorem 3, 



>(l + o(l))ifmes(' U He,k)n^"'e-"'/2. 

similarly, 

P(/,(T) > 6) < U ^^.k) > ft) 



kGC^ 

72 



<(l + o(l))ifmes(' U H,,k)^x^~'e-" 
Jordan measurability of T implies that 

mesf IJ He,kj-mesf J 'E^^A^^^. 

Therefore, 

P(la{T) > 6) = (1 + o(l))i/ mes(T)tx'^-ie-"'/^ □ 
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4. Proof for Theorem 2. In this section we present the proof of Theo- 
rem 2. We arrange all the lemmas and their proofs in the Appendix. 

Proof of Theorem 2. We evaluate the probability by conditioning 
on (/(0),9/(0),52/(0)), 



h{w,y,z) 



(12) 



P / e'"^Wdt>6 



f{0) = u-w,dfiO) = y, 

5^/(0) = ufiQ2 + z \ dw dy dz 



n 



h{w, y,z)p( / e^^W+'^^'W dt>b]dwdy dz, 



where 7^ = Ri<i+'^)i<i+2)/2 ^^^^ f^f^^^ jg density function of (/(O), 9/(0), 
5^/(0)) evaluated at {u — 'w,y,ufio2 + z). Now we take a closer look at the 
integrand inside the above integral. Conditional on /(O) =u — w,df{Q) = y, 

a2/(0) = Ufl02 + z, 



J\t 



^uE{t)+ug{t) 



|t|oo<£ 



/ exp< a 

■^|t|oo<e I 



1 



u-w + y'^t + -t^{-ul + z)t 



+ gz{t) + gA{t) + R{t) + g{t) 



dt 



■ det{ul — z) 



-1/2 



^a{u-w+{l/2)y'^ {ul~x)-^y} 



\{uI-7.)-^/H\^<e 



1 



X exp< a 



_{t-{uI-z)-^'^y)'{t-{uI-7.)-^/^y) 
+ g,({uI-zr'/H)+g^{{uI-z)-'/H) 

+ R{{ul - z)-^/h) + g{{ul - z)-'^/h) 



dt. 



For the second equality, we plugged in (5). For the last step, we first change 
the variable from t to {ul — z)^/^t and then write the exponent in a quadratic 
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form of t. We write the term inside the exponent without the factor a as 

Jit) = -Ut- i^i - ^r'^'yVit - i^i - ^r'^'y) + 93{{ui - zr'/H) 

(13) 

+ g^iiul - z)-^/h) + Riiul - zy'^h), 

which is asymptotically a quadratic form. But, as is shown later, and 
terms do play a role in the calculation. Also, it is useful to keep in mind 
that J{t) depends on y and z. Hence, we can write 

-'|i|oo<e 

= det(n/ - z)-i/2e-{«— +{i/2)yT(«/-z)-i.y} 



e 

l(n/-z)-i/2tl^<£ 



<7j(i)+ag({n/-z)-i/2j) 



Let 



(14) e^°=/ e-'^"'^^''''dt = {2^/a)'"\ 

JR'i 

Let u solve 

(15) ^-d/2gan+Ho ^ 5_ 



Then, 



if and only if 



dt > b, 



det{ul - z)-l/2g<T{n-«;+(l/2)yT(„/_z)-ij^} / ^aJit)+agi(uI-z)-^/H) 

We take the logarithm on both sides and rewrite the above inequality and 
have 

(16) 0< ^y'^ {ul - z)-^ y - aw - ^logdet{I - u~^z) 

+ log / e'^-^^*) dt - Ho 

i|(M/-z)-i/2i|^<e 
+ log^exp(cJ5((u/ - z)-^/^S)) 

(17) = Aiw,y, z) + logSexp(a<7((u/ - z)-^/^^)). 
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where S* is a random variable on the region that |('u/-z)~i/25|^ <£ ^itj^ 
density proportional to e'^'^^^^ and 



A{w,y,z) = ^y'^{ul -z) - aw - ^logdet{I - u ^z) 
(18) 

+ log / e'^-^^^Ut - Ho. 



'|(«7-z)-i/2j|<£ 

Thanks to Lemma 1, we only need to consider the set that 
£ = 11/(0) - n| <n2'5+^o, 

(19) 

Also, by abusing notation, we write 

(20) C = {\w\ < U^^+'\\y\ < < ^l/2+eo}_ 

Lemma 2 gives the form of h{w, y, z). We plug in the results in Lemmas 1 
and 2, 

p{Es) = / h{w,y,z) 

X P{A{w,y,z) +logEex.p{ag{{uI -zy'^^'^S)) > 0) dwdydz 
= o(l)'u""e-"'/2 

+ / h{w,y,z) 



ic 

X 



P{A{w,y,z) +\ogEexp{ag{{uI -zy^/^S))>0)dwdydz 



= o(l)u-°e 

(21) + - iri-^/2 



j P{A{w,y,z) + logEexp{ag{{uI - z)-^/^S)) > 0) 



X exp 



1 2 ^ N 

-u +-A{w,y,z) 

2 a 

+ ^y'^{I-(I-z/u)-')y 

"To"; ri r -Z fl22 ^ 

2 1- f-l20lJ^22 /^02 ^ 

+ ^ log det(/ — u~^z) 

--log/ e-^Wdt + -/7o 

0" yi(u/-z)-i/2tl<e cr 



dwdydz. 
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We define 



Z=\u^ + ^A{w, y, z) + \y^{I -{I- z/n)"i)y 

(22) + -\ ^ + -z^fi22 Z + logdet(/ - n'^z) 

--log/ e^'^Wdt + -i?o, 

y|(«/-z)-l/2t|<£ 0- 

and proceed with some tedious algebra to write X in a friendly form for 
integration. First notice that 

oo 
fe=0 

Plug this into the third term of I and obtain 

I=\u' + -A{w, y, z)-{l + 0{\z\/u))^y^zy 
la 2 

2 1 - ;U20^22 A*02 



-log/ e-^Wdt + !ii/o. 



J|(«/-z)-l/2t|<£ 

Situation 1 of Lemma 5. Adopt the notation in Lemmas 4 and 5. Note 
that according to the definition of Y in Lemma 4 that 

Y = {yf,i = 1, . . . 2yiyj, l<i<j< , 

d d{d-l)/2 

we obtain that 

y~^zy = Y~^ z. 

We plug in results of Lemmas 4 and 5. First, considering the first situation 
in Lemma 5, that is, £,i = jC n {| (u/ — z) ' y\oo 

we have 

I=lu^ + -Aiw,y,z) - (1 + Oi\z\/u))^Y^z 
la I 

^ z^l_,^^^z + — logdet{I-u ^z) 

2 1 - ;U20/^22 ^^02 ^ 

+ ^(u-'y + l/a)'^^i22iu-^Y + 1/a) - ^1^^*221 
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+ !^Ho - -H{{ul - z)-V2y, _ z)l/2^ ,) 
a a 

-^E^"-C'(0) + o(l). 

i 

Also, it is useful to keep in mind that 1 is NOT a vector of I's. The next 
step is to plug in the result of Lemma 3 and replace the logdet(/ — u~^z) 
term by 

-u~^ TV(z) + ^u-^£^{z) + o{u~^) 

and obtain 

X = 1 2 ^ !iA{uj,y,z) - (1 + 0{\z\/u))^Y^z 
Z a I 

2 1 - ^20/^22 A*02 2 2cj 4o-M 

+ l{u~^Y + l/a)^^22(^^~'^ + I/^t) - ^1^/^221 

+ - -Hiiul - zy^/^y, (ul - zf/\e) 
a a 



\Y.^niiC{Q)+o{l). 



8ct2 

Then, we group the terms —(1 + 0{\z\/u))^^Y^ z and —-^l^z and leave 
the 0{\z\/u) to the end and have 



1 2 . / N 1 + fJ'20f^22 



l^^2 



2 CJ 2 1- /i20M22 /^02 

+ ^^^^^22^^ - ^(^"^^ + 1/^)^^ 

+ liu-^Y + l/a)'^fi22{u-^Y + 1/a) 



+ ItHo - -H{{ul - z)-i/2y, (uI - zfl\e) 
a a 

-^2 l"'/^22l - A E ^--^(0) + 



8cj2 8a2 
+ 0(n-2|z|2|y|2) + 0(^-i£:2(2)). 
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Note that second and third hnes in the above display is in fact in a quadratic 
form. We then have 

la 2 1- /i20/^22 ^^02 



1 

+ 2 



+ !^Ho - -H{{ul - ^r'l'y, {ul - z)'/',e) 



-1/2 1 1/2. , -, I \ 

^22 ^ - 2^*22 ^ + 1/^^) 

-1/2 1 1/2. , 1 / ^ 
^22 o^22 



T 



a 



a 



-^2 l"'^*22l - ^ E ^--^(0) + 



8(j2 ^ 8a2 
+ 0(n-2|z|2|y|2) + o(^-i£:2(2)). 
Now, consider another change of variable, 
A = A{w,y,z), 



(23) 



y = y- 



Then, by noting that /i20 is a row vector in which the first d entries are — I's 
and the rest are O's, we have 



Therefore, we have 
2 a 2 



^ + /^20M22^^^-B + 7r-/W2ol + o(l). 
2a 



a 



+ 



1 {-A/a + mf^22^^B + (l/(2a))/i2ol + o(l))' 

2 1 — /X20^22^^02 
(T a 



We write X„ = Op(l) if — )• in probability as u — )■ 00. We insert the above 
form back to the integral in (21) and apply Lemma 7, 



p(H,)=o(l)n--e-«'/2 
1 

+ 



(27r)(d+l)('i+2)/4 



ir|-V2 
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X I P{u-A>Op{l)) 



X exp 



+ 



2 0- 2 

1 i-A/a + fi20^i22^^B + l/(2cj)/i2ol + 



1 — /^20A'22^/^02 

a 

^l^ml -^^diiiiC{0) + o(l) 



2 

u 
a 



8a^ 



+ 0{u-^\z\^\y\^) + Oiu-'£-'{z)) 



lc2/ 



I dw dy dz. 



Note that Jacobian determinant is 
det' ^(^'"'^^ 



'T-Met(/i22)^/2. 



Note that when \ (ul — z)~'^^'^y\ < ku^ — u^^"^ (the first situation in Lemma 5), 

uHq - uH{{uI - z)-^/2y, (mI - zf/"^ ,e) = o(l). 
Then, with another change of variable, A' = uA, the integration on Ci is 



ir|-V2 



(27r)('^+i)('^+2)/4 

X /" P(u-^>Op(l)) 



'£1 

X exp 



(24) 



+ 



2 o- 2 

1 (-A/(7 + fi20f^22^^B + 1/(2ct)m201 + o(l))' 



1 - /^20/^22^/^02 

+ ^i7,_!li7((^/_z)-l/2y,(^J_^)l/2^,) 
(T 



2 



8C72 



+ o(n-2|z|2|y|2) + 0(n-^£:^(z)) 



lc2/ 
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^-i|r|-i/2 



det(;U22)'/' 



(27r)('^+i)('^+2)/4 



X e 

X 



(l/2)n2+(l/(8ff2))iT^22l+(l/(8a2))j:^a,,,,C(0) 



/ Piu-A>Op{l)) 



'Ci 

X exp 



- ■ 

cj 2 



2(1 - Ai20/^22Vo2) 
+ 0(l) + 0(n-2|z|2|y|2+^,-l£:2(2)) 



+ 



dy 



-1/2 



■det(/Z22)'/'^x-i 



(27r)('^+l)('^+2)/4 
X g-(l/2)«2+(i/(8a2))lT^22l+(l/(8<x2)) 9™C(0) 



/ P{A'>Op{l)) 



X exp 



^' B^B 



A'/{au) + mt^22^B + l/(2a)/i2ol + o(l))^ 

2(1 - /i20^^2Vo2) 
+ 0(l)+0(n-2|z|2|y|2+^-l£:2(2)) 



+ 



d^' dB dy. 



The second equality is a change of variable from {w,y,z) to {A,B,y). The 
third equality is a change of variable from {A,B,y) to {A',B,y). Note that 
P(^' > Op(l)) /(A' > 0) as u — > oo. In addition, on the set C, 



0{u~^\z\'\y\' + u-'£\z)) = 0(7x-i+2'5+2-o|z|2) 



1/4 ^-l+25+2eoU|2 



By choosing 6 and small enough, when \B\ <u 
\B\ > n^/^, \B\ = e(|z|), therefore, 

^ + 0(n-l+25+2.o |^|2) ^ (1 + 0(1))^ 

The integrant in (24) has the following bound, for A' > 
P{A'>Op{l)) 



:o(l); 



X exp 



A' ^ B^B ^ {-A'/{au) + fi2of^^2^'^B + {l/{2a))fi2ol + o{l)f 

a 2 2(1 - /X20/i22Vo2) 
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<2exp|-i 



+ o{l) + 0{u-'\z\'\y\' + u-'8\z)) 



^ ^ ^ {-A'/{au) + l22ofi22^^B + (l/(2a))/i2ol)' 

a 2 2{1 - ^l20^J'22l^02) 



for 6' small enough. Note that the Op(l) is in fact — ulog-Eexp((7((Ml — 
z)~^/^5)). Thanks to the result of Lemma 7, the integral of the left-hand side 
of the above display in the region ^' < is o(l). By dominated convergence 
theorem, (24) equals 

^-i|r|-i/2 

(l+o(l))^27r)(rf+i)(rf+2)/4 

X det(/X22)^/^n"^e~(^/2)"'+(^/(^'"'))^''^22i+(i/(8^2))^^g^^^^C(0) 



(25) X / exp 

/a'>0,|j/|oo<km''+i/2-'u*/2+i/2 



A'/a + 



+ 



(/i20/i22'^'i3 + (l/(2cT))/i20l)' 



2(1 -/i20/^22V02) 

(1 + o(l))i? mes('uSe)'u-ie-(i/2V^ 



dA' dB dy 



where H is defined in (9). The above display is obtained by the fact that 
mes(nH£) = n'^mes(H£) = k'^u'^/'^^'^^ . 

Situations 2 and 3 of Lemma 5. For the second situation in Lemma 5, 
let £2 = Cn{Ku^ - < |(u/-z)-V2y|^ < (i ^j^g^.^ ^^.^^^^ 

ci > such that 



/ P{u-A>Op{l)) 

JC2 



X exp 



■2 ' -A + -B^B 



-u 



+ 



2 0- 2 

1 (-^/a + tl20l^22^^B + (l/(2cT))/i20l + o(l))' 



(26) 



2 1 - ^20^22^/^02 

+ - -i^((^/ - z)-i/2y, (n/ - z)V2, e) 

CJ 0" 

^l^fl22l - A 5^ 5iiiiC(0) + o(l) 



8a2 8cj2 

+ o{u-^\z\^\y\^) + 0{u-'£\z)) 
< {cief + o(l))i?mes(nH,)u-ie-(^/2)«2^ 
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For the third situation, C3 = Cn {{I + £i)ku^ < \{ul - zy^/^y\ < u'^+^"}, 
/ Piu-A>o^il)) 



X exp 



2 a 2 



+ 



1 {-A/ a + mi^22^^B + (l/(2a))/Z2ol + o{l)f 

2 1 - ^20A^22Vo2 



(27) 



+ !iHo - -H{{ul - ^r^r^y, {ul - z)^/^ e) 
a a 

- ^l^/i22l - J^5,,,,C(0) +0(1) 



+ o(n-^|zh2/n + 0(n-^£:^(z)) 



dy dz 



<0(1) 



^ -l^(l/2+5+eo)dg-(l/2)n2 



2, 

= o(l)mes(nH,)'u-ie^(i/2)«2_ 

We put (25), (26) and (27) together and conclude the proof. □ 

5. Proof for Theorem 3. Similar to Section 4, we arrange all the lemmas 
and their proofs in the Appendix. 

Proof of Theorem 3. Since the proofs for C"*" and C~ are complete 
analogue, we only provide the proof for C"*". We prove for the asymptotics by 
providing bounds from both sides. We first discuss the easy case: the lower 
bound. Note that 



>P(max7^(He,k) >6 
Vkec+ 

> PiU^eM) > ^) - E P(^-i^eM) > b,lA^e,k') > b). 

kGC+ kT^k' 



Thanks to Lemma 8, 



pU([j >&)>(! + 0(1)) Pil^^i^e, 



:)>&)• 
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I I I 

Fig. 2. Illustration of A, B and V. 

The rest of the proof is to estabhsh the asymptotic upper bound. To simpHfy 
our writing, we let 

(28) A = UE,), fi = /,(supH,,k), V = lJ [j H,,k'), 

where N is the set of neighbors of H^, that is, k G A/" if and only if 

inf |s-t|=0. 

An illustration of B and D is given in Figure 2. 
Further, let 

(29) bo = u~'-''/\ b-bo = {l- u-'-^/^)b = e-^^+'^W)--'-"'' 6, 
and uq solves 

and there exists cq > such that uo> u — cq logu. The first step in developing 
the upper bound is to use the following inequality 

P{A + B + V>b) 

< P{A >b-bo) + P{A <bo,A + B + V>b) 

(30) + P{bo<A<b-bo,A + B + V>b) 

< P{A >b-bo) + P{B + V>b-bo) 

+ P{A >bo,B + V>bo,A + B + V>b). 

From Theorem 2, 

P{A > 6 - 6o) = (1 + o{l))PiA > b). 
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The next step is to show that the last term in (30) is ignorable. Note that 
P{A >bo,B + V>bo,A + B + V>b) 

= P{A + B>b-bo,A>bo,B + V>bo,A + B + V>b) 
+ P{V >b-bo,A>bo,B + V>bo,A + B + V>b) 
+ P{A + B>bo,V>bo,A>bo,B + V>bo,A + B + V>b) 
<P{A + B>b-bo,A>bo,B + V>bo,A + B + V>b) 

+ 2P{V>bo,A>bo) 
= oiP{A>b)). 

The last step is due to Lemmas 9 and 10. By noting that #(C"*")ii~^~'^/2 = 
o{u~^), the conclusion of the theorem is immediate by induction, where #(•) 
is the cardinality of a set. □ 

APPENDIX: LEMMAS IN SECTIONS 4 AND 5 

Lemma 1 isolated the dominating event so that we will be in good shape 
to use Taylor's expansion. 

Lemma 1. There exist eo,5 > small enough and k large. Let e = 
Ku~^/'^^^ such that for any a > 0, 

p(^|/(0)-n| >'u2'5+^o or |9/(0)| >'u^/2+<5+eo 
|92/(0) - n/X2o| > u^/^+"\j H.e^W dt > b^ 

Proof. Note that there exists ci such that au < log6 + ci log log 6. Let 
au = log(6). Since we only consider the case that u is large, we always have 
mes(H£) < 1: 



p(^f{Q)<u-u^^^^'\ ef^^^dt>b 



< P(/(0) <u- 'u2'5+^o ^ sup /(t) > u) 

< CP{f{0) <u- u^^+''" \ sup f{t) > u)u'^-^e-^^/^. 

The last inequality is an application of Proposition 1. Because for any u' > u, 
for some ei > 0, 

inf e( f{t)\ sup f{t)=u'' 



.1^ „.2<5 



> u inf Cit) >u- K^eiu^\l + o(l)) 
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and 

sup Var(/(t)| sup/(i) = u') = 0{e^) = 0{u~^+^'), 

one can choose k large enough such that 

p(/(0) <u-u2^+^"|sup/(i) >u) =0(l)exp(-n^+^Vc2). 



Therefore, 



]u -°e-"'/2 



for all a > 0. Also, P(/(0) >u + u^^+^o) = o(l)ti-°e-(i/2)«2 ^ Hence, 
p[\f{Q)-u\>u^^+'\ ef^'Ut>h^ =o(l)n-°e-(l/2)«^ 
Similarly, we have 

p[\f{fd)-u\<u^^+'\\df{Q)\>u^''^+^+'\ e^'^'Utyh^ 

= o(l)n""e-(i/2K, 
p(^|/(0)-n| <n2'5+^M92/(0)-n/i2o| >'u^/^+'^y^ e-^W > 6^ 

The above two displays are immediate by noting that (/(O), 9/(0), 5^/(0)) 
is a multivariate Gaussian random vector. In addition, (/(O), 0^/(0)) is in- 
dependent of 9/(0) and the covariance between /(O) and 9^/(0) is /io2- D 

Lemma 2. Lei h{w,y,z) be the density o/ (/(O), 9/(0), 9^(0)). T/ien, 
/i('u;,y,z) 

^ |r|-i/2 



(27r)('^+i)('^+2)/4 ' 
X expl - + -^(w, y, z) + \y~^ {I -{I- z/u)-^)y 

2 1 - /X20/^22 /^02 ^ 



- log / e^-^(*) + nFo/cJ 

0" J\{uI-z)-^/^t\<e 
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1 - /^20M22^/^02 



/^22Vo2 



V 1 



^^20^^22 ^02 



_^_/f20^f22^_ 

-1 A^22Vo2/^20A^^2^ 

/-i22 1 

1 - ^^20/^22 ^02 / 



Proof. The form of T^^ in (31) is a result from linear algebra. The 
form of is direct application of the block matrix inverse from linear 
algebra. Note that 

^^^'^'"^= (2vr)('^+i)('^+^)/4 



1 



ir|-V2 



X exp^ --{u-w,z'^ + uij.2o,y 
By plugging in the form of T~^, we have 

{u - w)"^ 





y^y + 



1 - /^20/^22^/^02 
+ (z + Ufl02V 



-1 I /^22Vi02 ^20/^22^ 
f'22 "T 



2{u-w) 



/^20/^22 



1 — ^J'20^J■22 f^02 



1 - /^20^22^A'02 

(z + n;Uo2) 



(z + U/i02) 



« +y y — 2wu - 



-1 I /^22Vo2/^20Ai22^ 
A'22 ' 



1 — At20M22^/^02 

2; + 2tt; 



1 — /^20/^22^A'02 



A'20^22^ 



1 - A'20/^22^/^02 



2n 



+ — ^(w;,y,z)+y"^(/-(n7-z) ^)y + 



, (u; + /i20^22^ 



U^2 



(7 



1 — /^20/^22^/^02 



+ z'''//22^-2 H — logdet(/ — u ""^z) 



2n, 
— log 
a 



(u/-z)-i/2t|<e 

Therefore, we conclude the proof. □ 



a 
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Lemma 3. 

log(det(/ - u-^z)) = -u-^ Tr(z) + ^u-'^£^{z) + o{u-^), 

where Tr is the trace of a matrix, <f^(z) = Yli=i '^^^ ^''"^ eigen- 
values of z. 

Proof. The result is immediate by noting that 

d 

det(/ - u-^z) = Y[{1 - Xi/u) 

i=l 

and Tr(z) = Y^^^i Aj. □ 

Lemma 4. Let y = {yi, . . . ,yd)'^ and X ^ N{y/y/u,I/^/o■). Then, on the 
set C defined in (20), 

E{g;{X/^)+g^{X/^)) 
o 

i 

where 

Y = {yf,i = l,...,d, 2yiyj, 1 < i < j < d)^, 
1 = (1^_^,0^_^)T. 

d d{d-l)/2 

Proof. Using the derivatives in (4), we have that 

d 

dijkEiO) = -Y,dtjMC{0)m, 
1=1 

4fc£;(0) = {u + 0{\z\ + \w\))dijkiC{0). 
We plug this into the definition of g^ and 174 in (6) and obtain, on the set C, 
E{g3{X/V^)+g4X/V^)) 

= -l^~'^'Y.9fjkiC{0)E{X,X,X,yi) 

ijkl 

+ IT E dfj,iC{0)E{XajXkXi) + o{n~') 

ijkl 
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ijkl 

+ ^ E ^t^klCm(X^X,Xk{Xl - m/^)) + o{u~') 

ijkl 

E dtjkiC{0)y^y,ykyi - ^u~' ^ ym E 9unC{0) 

ijkl il j 

+ ^ E y^y^ E 4//^(o) + ^ E ^^'-^(o) + °(^^') 

jj / i 

-1^^^ E ^tjkiCiO)yiyjykyi - -^u"^ Yl y^y^ E ^tjiiC{0) 



4a 

ijkl ij 



8a2 

This last step is true because Y^iiyiyiJ2jdfijjC{0) = YlijyiyjYli9fjuC{0), 
which is just a change of index. Then, with the definition of Y and 1 in the 
statement of this lemma (note that 1 is NOT a vector of I's), we have 

E{g3{X/V^)+g4X/V^)) 

i 

= -^{U-'Y + l/cT)T/X22(n-^y + 1/(7) + ^1^221 

+ £iE^^"^(o) + °W- 

Lemma 5. Let J{t) be defined in (13). Then, on the set C the approxima- 
tions o/ i|(u/_z)-i/2i| e'^"'^*) (it under different situations are as follows: 

(1) W/ien |(u/-z)-l/2y|^ <Ky<5_y<5/2^ 

e'^-^^*) dt 



(«/-z)-l/2f|^<£ 

= exp[(T^((73(^/^/^) + 9a{X/V^)) 

where 

^H(y,^,e) ^ f ^-(a/1)(t-yV{t-y) 
J\T:'H\<e 
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and X is the random vector defined in Lemma 4- In addition, 
H{0, {ul - z)i/2^e) -Ho = o{u-^). 
(2) For any ei > 0, when ku^ - u^l'^ <\{ul - z)-^/2y|oo < (1 + £1)^"^ 



/ e'"'^(*) dt 

i|(w/-z)-l/2t|<£ 

< exp[(7£;(53(X/VI^) + g^{Xl^)) 

+ H{{ul - x)-^/^y, {ul - z.f'\e) + o{u-% 
(3) When (1 + ei)KU^ <\{ul - z)"^/2y|oo < 

/ e-^-^^*) dt 

J\{uI~z)-^/H\<e 

<^exp[aE{g3{X/V^)+g4{X/V^)) + Ho + o{u-')]. 

Proof. Note that 

f e'^'^^*) dt 

i|(u7-z)-i/2f|<e 

^ gH((M/-z)-i/2j/^(„/„z)-i/2^e) 

xE{eM^93{{uI-z)-^/^X') 

+ agiiiul - -lY^I^X') + aR[{ul - 'lY^I'^X% 

|(n/-z)-i/2x'| <e}. 
Also, (ul - •lY^I'^X' = (1 + 0{zlu))X' 1^ and 

X' = X- yj^ + {ul - 'zY^I'^y = X + 0(u-^l'^)\Ly\. 
For the first situation, \{ul — z)~^/^y|oo < k-u^ — u^^"^ and \z\ < u^/^+'^o, 
S{exp[cJ53((n/ - z)-i/2;^') + agi{{ul - zY^'^X') + cji?((n/ - zY^''^X'Y 

\{uI-zY^'^X'\<e} 



: E{exp[cJ53(X7V^) + ag^iX'/^) + aR{X'/^) + o( 



u 



-In 



|(u/-z)-i/2x'| <e} 
= £;{exp[a53(^7^) + (ygi{X'/^) + (ji?(X7VS) + o(m-1)]} 
= £;{exp[(T53(X/V^) + agYX/V^) + + o('u-i)]}. 

In addition, because \R{t)\ = 0{u^/'^+h^), then 

= S{exp[a(73(^/^/^) + CT54(^/V^) + o{u'^Y\- 
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Further, g'i{X/^) = o{u~^/'^+^) and gi{X/^/u) = o{u-^/^+^), then by re- 
peatedly using Talyor's expansion, we have 

E{exp[ag3{X'/^/^) + agi{X'/^) + aR{X' / ^)]-\{uI - 2.)-^'^X'\ < e} 

= eME[ag3{X/V^) + ag^iX/V^) + o{u-^)]}. 
For the second situation, the inequahty is immediate by noting that 
E{e^p[ag3{X'/^/^) + ag^iX'/^) + aR{X' /^/u)\-\{ul - 7.y^'^X'\ < e} 

< E{exp[ag3{X'/,/^) + ag^iX'/V^) + aR{X'/V^)]}. 

For the third situation, note that the integral is not focusing on the domi- 
nating part, and the conclusion follows immediately. □ 

The next lemma is known as the Borel-TIS lemma, which was proved 
independently by [16, 20]. 

Lemma 6 (Borel-TIS). Let f{t), t£U,U is a parameter set, be a mean 
zero Gaussian random field, f is almost surely bounded on U. Then, 

£'(sup/(t) ) < oo 
P( max/m - £;[max/(t)l > b) < e~^'/(H), 



and 



where 

f^w = maxVar[/(t)]. 

Lemma 7. Let log E ex.p{g{{ul — z)~^/'^ S)) be defined in (17), then there 
exists a A > such that for all x > 

P(n-3/2-35|logEexp(cj(7(('u/ - z)"^/^S))\ > x) < e"^^' 
for u sufficiently large. 

Proof. Note that g{t) is a mean zero Gaussian random field with 
Var(5(t)) = 0{\t\^) and |5| < ku^ . A direct application of the Borel-TIS 
lemma yields the result of this lemma. □ 



Lemma 8. For each k / k', 

P(/,(H,,k) > b, I^{E,y) >b) = 0(l)^x'='-ie-«'/2-QH. 
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Proof. Without loss of generality, we consider H^^k' = ^e- If and ^e,k' 
are connected to each other, 

P(/^(H,,k) > b,UEe) >b)< P(/,(H,UH,,k) > 2b) 

= 0(l)m(H,)n^-ie-(i/2)(.+iog2+o(i))2 

The second step is an application of Theorem 2. If and H^^k' are not 
connected, that is, inf^es^.teH^ \s — t\ > e = then 

< P( sup f{t) > u — clogu, sup /(t)>u — clogu 
<P( sup /(t) + /(s)>2u-2clogn 

" V ^4 - e(i)^i-i+2<5 y 

where Z is a standard Gaussian random variable. The last inequality is an 
application of the Borel-TIS lemma in Lemma 6. □ 

Lemma 9. Let A and T> he defined in (28) and 6o be defined in (29). 
Then, 

P{V >hQ,A> bo) = o{P{A > b)). 

Proof. Similar to the second case in the proof of Lemma 8, we have 
PiD>bo,A>bo) 

< p(^sup f{t) > u — cilogu, sup /(t) > M — ci logn 

^= UkGC+\{0,k} ^e.k 



<P\^ sup /(s) + > 2u- 2cilogti 

s6=e,tGUkgC+\{0,k} =e>k 
<^ag-(2«-2cilog«)V(2{4-2K2„-l+2*)) ^ g-«V2-e(l)«l+2* ^ 

The conclusion follows immediately. □ 

Lemma 10. Let A, B and V he defined in (28) and bo he defined in (29). 
Then, 

P{A + B>b-bo,A>bo,B + V>bo) = o{l)P{A > b). 
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Proof. Note that there exists c' > such that 
P{A + B>b-bo,A>bo,B + V>bo) 

< p(a + B > b - bo, sup f{t) >u-c'logu, 



sup f{t)>u — c'logu^. 

UkeC + \{0} ^e.k 

It suffices to show that the right-hand side of the above inequahty is 
o{l)P{A > b) and also o{l)P{A + B > b). In order to do so, we will bor- 
row part of the derivations in the proof of Theorem 2. Let solve 



b-ba. 



Note that because bo = u ^ '^/'^b, we have |m — u^^l = o{u and e "^/^ 
(1 + o(l))e-"*/2. By the results in (25), (26) and (27), we have 

p(^A + B > b — bo^svcp f{t) > u — d \ogu, sup f{t)>u — clogu 

UkGC+\{0} ^e.k 

= o{l)P{A + B>b-bo) 

+ (1 + o(l))a-^ det(r)-i/2 det(/i22)'/'tx-' 



X e 

X 



{|j/|oo<3Kni/2+«} 



P[A' > Op{l),supE{t) + g{t) > - c'logu*, 
sup E{t) + g{t) > — d log 

UkeC+\{0} 



X exp 



A! B^B 



{-J^l{au) + mi^22^^B + (l/(2a))/Z2ol 

2(1 -/i20/^^2Vo2) 



+ 



^2^ 



(iyl' d5 dy. 



Note that the only change in the above display from (25) is the probability 
inside the integral. In what follows, we show that it is almost always o(l). 
Note that Yar {g{t)) = 0{\t\^). Therefore, for any /(O) <u + u^" with eo < 
6/2, if 



(32) 



sup E{t) <u — clogu — Q{u ^) or 
sup E{t) <u — c'logu- Q{u^^), 
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then 

p(supE{t)+g{t) > u-clogu, sup E{t)+g{t) > u-c'loguj = o{l). 

UkeC+\{0} "e,k 

This fact imphes that g{t) can be basically ignored. Therefore, it is useful 
to keep in mind that "E{t) w /(t)." 
Since 

P(sup/(t) >u + M-^+"°) = o{l)P{A + B>b-bo), 

we only need to consider the case that supy < n + u~^~^'^°. Given the 
form 

E{t) = u^-w + y^t + \t^{-uj + z)t + g-i{t) + gi{t) + R{t), 

which is asymptotically quadratic in 'B^e, let t* = argsup^j^^ -^(^)- On the set 
that supT^ E{t) <u + u~^~^^° , we have 

sup E{t) < u — c logu — @{u^^). 

|i-i*|>2«-l/2 + e:o/2 

Let be the border of H^. Then 

sup E{t) > u — c' logu — Q{u~^) and 

(33) 

sup E{t) > u — c'logu — @{u 
only when init^dEe 1^ ~ < ^* This implies that 



inf |t— (li*/ — z) < u. 



-1/2+eo 



Therefore, t* = argsup/(t) must be very closed to the boundary of so as 
to have (33) hold. 

Therefore, for all Eq < S 

p(A + B>b- 5o,sup/(t) >u- c'logu, 



sup E{t) + g{t) > u — c'logu 

UkeC+\{0} ^e.k 

= o{l)P{A + B>b-bo) 
+ (1 + o(l))a"i det(r)-i/2 det(/i22)'/'n-i 

X / P(^'>0p(l), 
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sup-E(t) + g{t) > u — c'logu, 



sup E{t) + g{t) > u — c' log u 
Ukec+\{o} "E.k 



X exp 



A' B^B 



{-A'/{au) + l22of^22^^B + (l/(2a))//2ol) 

2(1 - /X20/^^2Vo2) 

= o(l)P(^ + e>&-6o)- 
The last equation is because 

mes^jy :^inf \t — {u^I — z)"^y| < u^^^^"*"^"!^ 

= o(mes({y -Mul- z)-i/2y|^ < _ ^&/2^^y 

Hereby, we conclude the proof. □ 



+ 



dA'dBdy 



[1 
[2: 

[3 

[6: 

[7] 
[9 

[lo; 
[11 



REFERENCES 

Adler, R. J. (1981). The Geometry of Random Fields. Wiley, Chichester. 
MR0611857 

Adler, R. J., Blanchet, J. H. and Liu, J. C. (2008). Efficient simulation for 

tail probabilities of Gaussian random fields. In Proceedings of the 40th Winter 

Simulation Conference {Miami, Florida) 328-336. 
Adler, R. J., Blanchet, J. H. and Liu, J. C. (2009). Efficient Monte Carlo for 

large excursions of Gaussian random fields. Preprint. 
Adler, R. J., Samorodnitsky, G. and Taylor, J. E. (2009). High level excursion 

set geometry for non-Gaussian infinitely divisible random fields. Preprint. 
Adler, R. J. and Taylor, J. E. (2007). Random Fields and Geometry. Springer, 

New York. MR2319516 
Adler, R. J., Taylor, J. E. and Worsley, K. J. (2009). Apphcations of RANDOM 

FIELDS AND GEOMETRY: Foundations and case studies. Preprint. 
Ahsan, S. M. (1978). Portfolio selection in a lognormal securities market. Journal 

of Economics 38 105-118. 
Asmussen, S. and Rojas-Nandayapa, L. (2008). Asymptotics of sums of lognormal 

random variables with Gaussian copula. Statist. Probab. Lett. 78 2709-2714. 

MR2465111 

AzAlS, J.-M. and Wschebor, M. (2005). On the distribution of the maximum of 
a Gaussian field with d parameters. Ann. Appl. Probab. 15 254-278. MR2115043 

AzAlS, J.-M. and Wschebor, M. (2008). A general expression for the distribution of 
the maximum of a Gaussian field and the approximation of the tail. Stochastic 
Process. Appl. 118 1190-1218. MR2428714 

AzAlS, J.-M. and Wschebor, M. (2009). Level Sets and Extrema of Random Pro- 
cesses and Fields. Wiley, Hoboken, NJ. MR2478201 



TAIL APPROXIMATIONS OF INTEGRALS OF GAUSSIAN RANDOM FIELDS 35 



[12] Basak, S. and Shapiro, A. (2001). Value-at-risk-based risk management: Optimal 
policies and asset prices. Review of Financial Studies 14 371-405. 

[13] Berman, S. M. (1985). An asymptotic formula for the distribution of the maximum 
of a Gaussian process with stationary increments. J. Appl. Probab. 22 454-460. 
MR0789369 

[14] Black, F. and Scholes, M. (1973). Pricing of options and corporate habilities. 
Journal of Political Economy 81 637-654. 

[15] Blanchet, J., JuNEJA, J. and Rojas-Nandayapa, L. (2012). Efficient tail estima- 
tion for sums of correlated lognormals. Ann. Oper. Res. To appear. 

[16] BORELL, C. (1975). The Brunn-Minkowski inequality in Gauss space. Invent. Math. 
30 207-216. MR0399402 

[17] BORELL, C. (2003). The Ehrhard inequality. C. R. Math. Acad. Sci. Pans 337 663- 
666. MR2030108 

[18] Campbell, M. J. (1994). Time-series regression for counts — an investigation into 
the relationship between sudden-infant-death-syndrome and environmental- 
temperature. J. Roy. Statist. Sac. Ser. A 157 191-208. 

[19] Chan, K. S. and Ledolter, J. (1995). Monte Carlo EM estimation for time series 
models involving counts. J. Amer. Statist. Assoc. 90 242-252. MR1325132 

[20] Cirel'son, B. S., Ibragimov, I. A. and Sudakov, V. N. (1976). Norms of Gaussian 
sample functions. In Proceedings of the Third Japan-USSR Symposium on Prob- 
ability Theory (Tashkent, 1975). Lecture Notes in Math. 550 20-41. Springer, 
Berlin. MR0458556 

[21] Cox, D. R. (1955). Some statistical methods connected with series of events. J. Roy. 

Statist. Soc. Ser. B. 17 129-157; discussion, 157-164. MR0092301 
[22] Cox, D. R. and Isham, V. (1980). Point Processes. Chapman and HaU, London. 

MR0598033 

[23] Davis, R. A., Dunsmuir, W. T. M. and Wang, Y. (2000). On autocorrelation in 
a Poisson regression model. Biometrika 87 491-505. MR1789805 

[24] Denisov, D., Dieker, A. B. and Shneer, V. (2008). Large deviations for random 
walks under subexponentiality: The big-jump domain. Ann. Probab. 36 1946- 
1991. MR2440928 

[25] Deutsch, H.-P. (2004). Derivatives and Internal Models. Finance and Capital Mar- 
kets, 3rd ed. Palgrave Macmillan, Hampshire. 

[26] Duffie, D. and Pan, J. (1997). An overview of value at risk. The Journal of Deriva- 
tives 4 7-49. 

[27] Dufresne, D. (2001). The integral of geometric Brownian motion. Adv. in Appl. 
Probab. 33 223-241. MR1825324 

[28] FOSS, S. and Richard, A. (2010). On sums of conditionally independent subexpo- 
nential random variables. Math. Oper. Res. 35 102-119. MR2676758 

[29] Glasserman, p. (2004). Monte Carlo Methods in Financial Engineering. Applica- 
tions of Mathematics (New York) 53. Springer, New York. MR1999614 

[30] Glasserman, P., Heidelberger, P. and Shahabuddin, P. (2000). Variance re- 
duction techniques for estimating value-at-risk. Manag. Sci. 46 1349-1364. 

[31] Landau, H. J. and Shepp, L. A. (1970). On the supremum of a Gaussian process. 
Sankhyd Ser. A 32 369-378. MR0286167 

[32] Ledoux, M. and Talagrand, M. (1991). Probability in Banach Spaces: Isoperimetry 
and Processes. Ergebnisse der Mathematik und Ihrer Crenzgebiete (3) [Results 
in Mathematics and Related Areas (3)] 23. Springer, Berlin. MR1102015 

[33] Marcus, M. B. and Shepp, L. A. (1970). Continuity of Gaussian processes. Trans. 
Amer. Math. Soc. 151 377-391. MR0264749 



36 



J. LIU 



[34] PiTERBARG, V. I. (1996). AsymptoUc Methods in the Theory of Gaussian Processes 

and Fields. Translations of Mathematical Monographs 148. Amer. Math. Soc, 

Providence, RI. MR1361884 
[35] SUDAKOV, V. N. and Tsirelson, B. S. (1974). Extremal properties of half spaces 

for splierically invariant measures. Zap. Nauchn. Sem. LOMI 45 75-82. 
[36] Sun, J. (1993). Tail probabilities of the maxima of Gaussian random fields. Ann. 

Probab. 21 34-71. MR1207215 
[37] Talagrand, M. (1996). Majorizing measures: The generic chaining. Ann. Probab. 

24 1049-1103. MR1411488 
[38] Taylor, J., Takemura, A. and Adler, R. J. (2005). Validity of the expected Euler 

characteristic heuristic. Ann. Probab. 33 1362-1396. MR2150192 
[39] YOR, M. (1992). On some exponential functionals of Brownian motion. Adv. m Appl. 

Probab. 24 509-531. MR1174378 
[40] Zeger, S. L. (1988). A regression model for time series of counts. Biometrika 75 

621-629. MR0995107 



Columbia University 
1255 Amsterdam Ave. 
New York, New York 10027 
USA 

E-MAIL: jcliu@stat.coIumbia.edu 



